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In a hadronic gas with three conserved charges (electric charge, baryon number, and strangeness) 
we employ the hadron resonance gas model to compute both diagonal and off-diagonal susceptibili¬ 
ties. We model the effect of chemical freeze-out in two ways: one in which all particle numbers are 
conserved below the chemical freeze-out temperature and one which takes into account resonance 
decays. We then briefly discuss possible implications these results may have on two active areas of 
research, hydrodynamic fluctuations and the search for the QCD critical point. 


I. INTRODUCTION 

Analysis of fluctuations has been of great interest 
to the heavy ion collision physics community in recent 
years. One line of research examines the impact of pre¬ 
equilibrium initial state fluctuations on observed final 
state momentum anisotropies (see [1] for a review). In 
addition to initial state effects, fluctuations of a thermal 
nature will occur after the system formed in a heavy ion 
collision has attained local equilibrium. The enhance¬ 
ment of such fluctuations is an expected experimental 
signature of the QCD critical point [2-4]. The principle 
goal of the Beam Energy Scan at RHIC is the discovery 
of the QCD critical point, and the STAR collaboration 
has recently reported results in this direction [5]. 

Thermal fluctuations can also lead to observable par¬ 
ticle correlations [6-8]. These correlations can be ex¬ 
tended over a large pseudorapidity range because they 
are sourced by fluctuations which propagate or diffuse 
according to the equations of hydrodynamics. Such cor¬ 
relations are experimentally studied by measuring bal¬ 
ance functions [9, 10]. 

To examine the effect of (thermal) conserved charge 
fluctuations over the history of a heavy-ion collision, one 
needs thermodynamical input, as the magnitude of such 
fluctuations are controlled by the susceptibilities 
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Xa/3 = o 7 , (1) 
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where P is the pressure, /i denotes chemical potential, 
and the abstract indices a , (3 denote the particular con¬ 
served charge under consideration. In this work, a can 
take on the value Q, B, S for electric charge, baryon num¬ 
ber, or strangeness. (In fact, as explained in [8], the rel¬ 
evant quantity is actually XapT/s where s is the entropy 
density.) One obvious source of input on these thermo¬ 
dynamic quantities is lattice QCD, as the susceptibilities 
have been calculated there [11-13]. However, as the lat¬ 
tice assumes full thermal and chemical equilibrium, its 
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results may not accurately reflect the thermodynamics 
in a heavy ion collision where chemical equilibrium is not 
always maintained. It is the effect of chemical freeze-out 
on the conserved charge susceptibilities which we wish to 
explore in this work. For the highest energy collisions 
at RHIC and LHC, the net charge, baryon number, and 
strangeness of the fireball are approximately zero. Hence, 
for the remainder of this paper we calculate susceptibili¬ 
ties at fj-B = UQ — Us = 0. 

The effect of chemical freeze-out on thermodynamics in 
heavy ion collisions has a long history going back more 
than two decades [14-18]. However, these discussions 
have exclusively focused on pressure P and energy den¬ 
sity e, with emphasis the equation of state P(e) as this 
is necessary input for hydrodynamic simulations. It was 
found that despite the fact that both P(T ) and e(T) 
are both modified after including the effects of chemical 
freeze-out, P(e) is affected only slightly [16]. 

Given the recent interest in fluctuations we thus wish 
to reopen the line of inquiry with particular attention 
paid to the susceptibilities. To accomplish this, we em¬ 
ploy the hadron resonance gas model (HRG) [19-21] 
which (in full equilibrium) has proved to be an excellent 
approximation to lattice calculations at low temperatures 
(T < 180 MeV). By including the effects of chemical non¬ 
equilibration in the HRG, we expect the results to more 
accurately represent the thermodynamics of the system 
formed in a heavy ion collision. 

The following is the outline of the paper. In Sec. II, 
we review the thermodynamics of the HRG model. We 
then proceed to detail two ways of modeling chemical 
freeze-out. We refer to the first as “Full chemical Freeze 
Out” (FFO) and explain the details in Sec. Ill wherein 
all number changing processes cease below the freeze-out 
temperature T < T c h and hence all particle numbers are 
constant. We also provide some new analytical formulas 
for the chemical potentials and susceptibilities. In Sec. 
IV, we detail the second model of chemical freeze-out 
referred to as Partial Chemical Equilibrium (PCE) [14] 
which permits resonances with short lifetimes to decay. 
We present the results for the susceptibilities in Sec. V. 
We comment on interesting features of the results, possi¬ 
ble phenomenological implications and directions for fu¬ 
ture work in Sec. VI. 
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II. REVIEW OF HADRON RESONANCE GAS 
THERMODYNAMICS 

At temperatures below the deconfinement temperature 
(T < 180 MeV), the thermalized matter left in the wake 
of a heavy ion collision can be modeled as a gas of non¬ 
interacting, point-like hadrons. The effect of interactions 
is incorporated by including resonances. This model is 
parameter-free, and does a remarkably good job of ap¬ 
proximating lattice QCD data at sufficiently low tem¬ 
peratures [13]. 


Note that the total pressure is a function of the set of all 
chemical potentials, which we denote {/r}. The HRG is 
parameter free, but one must choose which particles to 
include in a given calculation. We include particles and 
resonances fisted by the particle data group (PDG) [22] 
with masses less than or equal to 2 GeV. More details on 
our included particles and their decays can be found in 
Appendix A. 


B. Conserved Charges and Susceptibilities in Full 
Equilibrium (FE) 


A. Pressure, Entropy Density, Number Density 


To access relevant thermodynamic quantities, one may 
start by considering the partial pressure of a given parti¬ 
cle species (labeled with the subscript ‘i’). Assuming the 
momentum distribution is isotropic, the expression is 
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where the upper/lower signs refer to fermions/bosons, 
E 2 = p 2 + m 2 and g.i is a spin degeneracy factor, We 
have included a chemical potential Hi associated with 
this particular particle. It is more convenient to write 
Hi in terms of chemical potentials associated with the 
conserved charges: baryon number (B), electric charge 
(Q), and strangeness (S) as we will see shortly. 

The number density of this particle is found by 


(T, Hi) 
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The second equality is a more familiar expression for the 
number density, and can be found by first taking the 
derivative of (2) and then performing an integration by 
parts. 

We will also need an expression for the partial entropy 
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The second equality can be found by first taking the 
derivative of (2) and then performing integration by 
parts, or more directly by using the thermodynamic iden¬ 
tity 


TSi — £i T Pi Hi^'i (b) 

where Si is the energy density of the particle “i”. In order 
to compute the total entropy density s and/or the total 
pressure P, one must sum over all hadronic species i , 

P(T,{/x}) = ^P i (T,/x i ). (6) 


We are interested in fluctuations of conserved charges; 
as such it pays to work with chemical potentials asso¬ 
ciated with these conserved charges only. In later sec¬ 
tions, when we consider chemical freeze-out, certain par¬ 
ticle numbers are conserved and hence we will introduce 
additional conserved “charges” and chemical potentials. 
The total density of charge a in the system is 

n a (T, {/i}) = TaM^Hi) (7) 

i 

where a,; is the conserved charge a of the Ath particle. 
For example, the electric charge density can be found by 
setting the abstract index a = Q: 

n Q(T, {h}) = ^ Q l n i (T, Hi) (8) 


Every conserved charge has a corresponding chemical 
potential, so by definition we could alternatively write 



Comparing (7) with (9) we see that 


( 9 ) 
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This can only be satisfied if the “particle” chemical po¬ 
tentials are related to the conserved charge chemical po¬ 
tentials as 


Hi / ) (%iHot 

OL 

or, in less abstract notation, 


( 11 ) 


Ah — Qid’Q T BiHB T SiHs (12) 

where Qi,Bi,Si are the electric charge, baryon number, 
and strangeness of the ith particle. 

We now consider the (3 x 3) matrix of susceptibilities 

dn a {T,{H}) o dn i{T,Hi) 


(13) 
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In terms of thermodynamic integrals, the components of 
the susceptibility matrix are 


X a p{T,{p}) 


E giOtifii 

2tt 2 T 
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r°° e( Ei fJ,i ' ) / T p 2 dp 
Jo [ e (E i -^)/T ± I] 2 


(14) 


For example, the baryon-strangeness susceptibility is 
found by setting a = B, /3 = S, 


Xbs (T,{p}) = Yl 


Qi Bi Si 


e {Ei ^i)/Tp2^p 


2 k 2 T J o [ e ( Ei -^)/T ±l y 


(15) 


C. Boltzmann Approximation 


If quantum statistics are necessary, one needs to use 
the expressions given in the previous subsection to com¬ 
pute relevant thermodynamic quantities. However, for 
physical hadronic masses and temperatures below decon¬ 
finement, quantum statistics represent only a small cor¬ 
rection to the classical results. The classical expressions 
are advantageous due to their analytical tractability. We 
will exclusively use the Boltzmann approximation for the 
remainder of this work. For Boltzmann statistics, we can 
neglect the ±1 in the distribution function. The remain¬ 
ing integrals can be done analytically leading to 


= nf E (T)e Mi/T 


(T, /J,i ) — Tli (T, Pi ) 
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T K 2 (rm/T) 


tk 
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where K 2 {x) is a modified Bessel function. The function 
n EE (T) is the number density in full equilibrium (FE), 
and where all chemical potentials vanish ps = Pq = 
ps = 0 (which is the case shortly after thermalization in 
a high energy heavy ion collision). 



The total entropy in full equilibrium (and where all chem¬ 
ical potentials vanish) is 


FE 


(T)=5> ee (T) 


mrii K 3 (rrii/T) 


T Kzim^T) 


The susceptibilities are given by 
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Given a list of hadrons with masses m,;, the (full equi¬ 
librium) susceptibilities as a function of temperature can 
be calculated using (20). Calculations of this sort were 
carried out and compared with lattice QCD results in 

[ 13 ]- 


III. FULL CHEMICAL FREEZE OUT (FFO) 

The results of the previous section are applicable for 
a system in full chemical equilibrium. The system cre¬ 
ated in a heavy ion collision is not static and has a finite 
size; it does not remain in chemical equilibrium through¬ 
out its entire evolution. Eventually, number changing 
processes freeze out, and the only remaining interactions 
between hadrons are elastic scattering processes. During 
this phase, the system is chemically frozen out (but still 
in local thermal equilibrium). 

A first attempt to model this phase can be made by 
considering the total number of each species of hadron 
(TVj = riiV where V is the volume of the system) to 
be fixed [15]. We use the superscript “FFO” to denote 
“Full Freeze-Out”, meaning each hadron species has a 
fixed number after chemical freeze-out. In reality, some 
particles decay; we include this effect in the next section. 

To maintain this chemical freeze-out condition, we 
must introduce additional chemical potentials corre¬ 
sponding to the new conserved particle numbers. As¬ 
suming each species of hadron is conserved after chemi¬ 
cal freeze-out, there will be one new chemical potential 
for each hadron. Hence: 


— QihQ + Bips + Sips + pf F ° (21) 

The additional chemical potential p FFO vanishes while 
the system is in chemical equilibrium, but becomes active 
after chemical freeze-out. One can think of this as sud¬ 
denly associating a conserved “charge” with each species 
of hadron. For example, 7r + particles have a (+1) 7r + 
“charge”, 7r° particles have a (+1) 7r° “charge” (and zero 
7r + ,7r - “charge”), etc... 

In order to remove any dependence on the volume of 
the system, it is more convenient to impose the condition 

ta(T, mF F °) = nJ E (T ch ) 

s(T,{p FE O}) s FE (T ch ) 1 J 

where the right hand side is a constant. The temper¬ 
ature T c h is the chemical freeze-out temperature below 
which number changing processes are no longer effective. 
This condition is equivalent to W = constant because 
(neglecting dissipative corrections) the total entropy sV 
is also conserved. 


A. Analytical Results 

The central problem is how to determine the temper¬ 
ature dependence of all of the chemical potentials pJ FO 
such that condition (22) is satisfied. In this subsection, 
we give some analytical formulas for the chemical poten¬ 
tials and susceptibilities which (to the best of our knowl¬ 
edge) have not previously appeared in the literature. 

Assuming ps = Hq = ps = 0 we have pj = p J FO - If 
we have a system of N to t species of particles, there are 
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Ntot ~ 1 constraints of the form 


B. Approximate Chemical Potentials for m; T 


n j(T, mJ fo ) nJ E (r ch ) 
7 l(EmF fo ) «, FE (T ch ) 


Substituting (16) and (18), one can solve for the differ¬ 
ence 


/iJ FO (T) - p FFO (T) f K 2 (m i /T)K 2 (m j /T cb y 
T n K 2 {m i /T ch )K 2 {m j /T) 

(24) 

Note that at T = T c h, the difference vanishes implying all 
of the chemical potentials are the same at that point (of 
course, by definition they should all vanish at T = T c h 
which we will show shortly). 

The final constraint is found by summing (17) over all 
hadrons and dividing both sides by s, which leads to 

n,(7>f°) 

^s{T,{p FF ° } ) 


rrij K 3 (mj/T) 

~T 1 ( 2(171 j/T) T 


u FFO 

n 


(25) 


In our calculations, we use (26) to compute the chem¬ 
ical potentials in the case of full freeze-out. However, to 
develop our intuition, it is useful to consider an approxi¬ 
mation which leads to simpler results. The HRG model 
should be applicable for temperatures of the order 100 
MeV < T < 180 MeV, below deconfinement and above 
kinetic freeze-out. For these temperatures, most hadrons 
have ml$> T except for the lowest mass ones (pions). 

For m T, one can use the asymptotic approximation 
for the modified Bessel functions 


K n (x) « 



1 + 0 



(31) 


If we assume that rrii T c3l > T , and keep only the 
leading order in the asymptotic expansion, the chemical 
potentials are found to be approximately linear in to with 
a temperature dependent offset 


Substituting (22), the expressions for rij and s in full 
equilibrium (18), (19), and the constraint (24) one can 
solve analytically for p FFO 


Mf F °(T) 

T 


= In 


A' 2 (TOi/T ch ) 


L K 2 (m i /T) J 


- A(T) (26) 


where A (T) is a temperature dependent offset which is 
independent of the mass to,; 


A (T) = 


and we have defined 


E j 9j K 2 {rrij /T c h) 


g»)-g(S) 


E j gjm/jK 2 (rn^/Tch) 


G{x) = Nfyy + ln i K ^ x )\ 

K2{X) 


(27) 

(28) 


Plots of the chemical potentials as a function of T are 
shown in Fig. 1 assuming T c h = 160 MeV. 

Once T c h is specified, and a set of particles to include in 
the HRG is fixed, one can use (26) to determine the chem¬ 
ical potential for each species of particle as a function of 
T. When these are known, one can go back and compute 
the susceptibilities by inserting the correct /.t FFO (T) for 
each species of particle, 

xl F p°(T) = M FFO (T)) (29) 


which can be written 


X^ U (T) 


T 2 e -A(T) 

2w 2 



(30) 

This analytical formula for susceptibility (and associated 
formulae for the chemical potentials) in “full freeze-out” 
is one of our main results. 


in (FE) +A(T) ■ 

(32) 

A similar approximation was used in previous calcula¬ 
tions [16]. Unfortunately, a simple approximation for A 
is more difficult to come by since the sum over j neces¬ 
sary to compute A runs over all hadrons, some of which 
have m ~ T. 

Assuming a freeze-out temperature of T c h = 160 MeV, 
for high mass particles (to ~ 1 GeV) and low temper¬ 
atures (T ~ 100 MeV), the approximate formula differs 
from the analytic result by less than 3%. 

In summary, for physical hadron masses and temper¬ 
atures realized in a heavy ion collision, the “full freeze 
out” chemical potentials are approximately linear with 
temperature (plus corrections). Corrections to this ap¬ 
proximation are due to both the contributions at to ~ T, 
and quantum statistics. 


A 


FFO 


(T) 


1 - 


T 

t7 


ch 


-T 


IV. INCLUDING DECAYS - PARTIAL 
CHEMICAL EQUILIBRIUM (PCE) 

In reality, the number of each particle is not precisely 
conserved after chemical freeze-out due to the fact that 
resonances decay. For example, the number of p mesons 
is not constant, since they will primarily decay into pi¬ 
ons. In order to account for this, we use the model of 
partial chemical equilibrium detailed in [14] which we re¬ 
view here for completeness. 

We consider only a subset of particles to be stable 
(i.e. those with lifetimes much longer than timescale of a 
heavy ion collision). We denote number of stable parti¬ 
cles and the total number of particles as Ns, N to t respec¬ 
tively. We follow [16] in choosing the stable hadrons to 
be 7 r, K, 77 , ui, N, r {, (j), A, E, S, fl (we include all members 
of any isospin multiplet and all antibaryons). Each of 
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FIG. 1. (Color Online) Results for the particle chemical potentials pf FO and pf CE which enforce the chemical freeze-out 
condition. Shown are the chemical potentials for the four lightest mesons (a) and the four lightest baryons (b). The dashed 
lines are the chemical potentials given the FFO freeze-out condition, the solid lines are the chemical potentials given the PCE 
freeze-out condition. We assumed T c h = 160 MeV. 


these stable particles has an associated chemical poten¬ 
tial pf GE - The generalization of (12) is now 

^r CE = Qi^Q + Bi^B + Sins + ^2 d ijP % CE -( 33 ) 

j, stable 

Note that the summation runs only over stable parti¬ 
cles. The coefficients dij are the “charges” associated 
with the stable conserved chemical potentials. For ex¬ 
ample, whereas in the full freeze-out case, a p + particle 
would have (+1) “ p + charge”, in the PCE case a p + par¬ 
ticle has (+1) “7r + charge” and (+1) “7r° charge” since a 
p + decays to 7T+7T 0 100% of the time. 

More generally, the d l:] coefficients are found by ex¬ 
amining the decay rates of particle i which result in one 
(or more) stable particles j. Computation of the dij re¬ 
quires branching ratios for each hadron and resonance 
included in the HRG. In Appendix A we provide more 
details on the determination of these decay rates, many 
of which are not measured precisely. The precise value 
d^ is found by multiplying the branching fraction for the 
decay by the number of stable particles j formed in that 
decay, and finally summing over all decay modes of i. In 
other words, dij is the average number of stable particles 
j formed from a decay of particle i. A stable particle can 
be thought of one which decays into itself 100% of the 
time, so for example d T o, = S^Oi- 

It is helpful to envision a matrix, d, which has N rows, 
and Ng columns. In appendix B we provide a portion of 
this matrix, (and hence give a subset of the dij coeffi¬ 
cients). 

The quantity which is now conserved after chemical 
freeze-out is 

Ni = Y J Njd ji (34) 

j 

which is the “effective number of stable particles”, in¬ 
cluding those hidden inside the unstable particles which 


will eventually decay. For example, if we had a system 
of only p and 7r particles, the number 


A%+ — N-+ + N„+ 


(35) 


would be conserved after chemical freeze-out, since a p + 
effectively counts as one tt + . Similarly, the quantity 


N^o — N^o + N p + + N p - 


(36) 


is conserved since both p + and p~ effectively count as 
one 7r°. 

As before, to remove any dependence on the volume 
of the system, the freeze-out condition is implemented 
using intensive quantities 


ni{T, {/I PCE }) 
S (T,{ A l p CE }) 


-FE 

n. 


C FE 


CTch) 

(Tch) 


= a 


(37) 


This criteria is more complicated than the full freeze-out 
case, since n l now depends on multiple chemical poten¬ 
tials. As before, we set up a system of Ns algebraic equa¬ 
tions for {p PCE }. There are Ns — 1 constraints which can 
be written without any reference to s: 


™i{T, {p PCE }) 

rjPCE 


C, 


nn(T,H PCE ) C n 


(38) 


It is advantageous to single out the heaviest stable parti¬ 
cle (Q) because it is both stable (dm = dm), and none of 
the other particles we consider decay into it (din = din). 
Hence, fin(T, {p PCE }) = nn(T, pn CE ). To make the de¬ 
pendence on the chemical potentials explicit, we can re¬ 
write these equations as 


C l nf E (T)e 


-PCE 

= Cn X! n J E ( T ) d A ex P 

j 


J , .-.PCE 
a jkP k 

. k 


T 


( 39 ) 

























6 


(a) 



T [MeV] 
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FIG. 2. (color online) Dimensionless thermodynamic quantities pressure (a), energy density (b), and entropy density (c) in 
the case of full chemical equilibrium (FE), full freeze-out (FFO), and partial chemical equilibrium (PCE) assuming T c h = 160 
MeV. The red points are lattice data adapted from [11] 


One more equation is necessary to close the set. We take 
it to be 


* s(T,{q PCE }) 
~ Q nn(T, p,Q GE ) 


(40) 


Substituting (17), and making use of (38), we find 


E 


s EE (T) 

FE frn\ eX P 


n 


'CO 


E 

_ k 


dikvl CE 

T 


„PCE 

>‘n 

= e T 


1+ E 


Citi 


PCE 


T 

(41) 


iVtot 

J2rf CE m= E & 

i —1 i, stable 


'“ PCE fL. 


The right hand side of (41) makes use of the fact that 

(42) 


The second sum runs only over the stable particles. This 
relation is straightforward to show given (33) and (34). 

Eqs. (39) and (41) form a nonlinear system of Ng alge¬ 
braic equations for the set of chemical potentials {/i PCE }. 
Unlike the “Full Freeze-Out” case, we have been unable 
to solve this system analytically. The system can be 
solved numerically (for example) with a matrix formu¬ 
lation of Newton’s method, but in practice we find the 
most efficient way to solve the system is with Mathe- 
matica’s built-in FindRoot [ ] function [23]. Results for 
both light mesons and baryons are shown in Fig. 1 


V. RESULTS AND DISCUSSION 

With the results for the chemical potentials of each 
particle in hand, we can then compute thermodynamic 
quantities of interest. The energy density, pressure and 
entropy density are shown in Fig. 2. We also include 
lattice data in our plots; we re-emphasize that the PCE 
and FFO results are not expected to agree with the lat¬ 
tice, since lattice calculations presuppose full chemical 
equilibrium. 

The susceptibilities can be seen in Fig 3. Note that 
the inclusion of chemical freeze-out can have a large ef¬ 
fect compared to a fully equilibrated system, especially 


at low temperatures where the full equilibrium suscep¬ 
tibilities are approximately vanishing. With regard to 
the inclusion of decays via PCE, we see that with only 
a few exceptions, the PCE scenario amounts to roughly 
an 0(10%) correction to the FFO scenario, as long as 
T > 100 MeV. 

When considering the magnitude of hydrodynamic 
fluctuations, the quantity \T/s is the relevant one [8]. 
We plot this combination in Fig. 4. At temperatures 
below T c h, Xa /3 °T/s is exactly constant (since x ~ Yli 
and Ui/s is constant). The inclusion of partial chemical 
equilibrium amounts to a very small correction in most 
cases. 

The baryon-baryon susceptibility is especially interest¬ 
ing. In the case of Xbb/s , the PCE and FFO curves 
lie on top of each other, and (within our numerical pre¬ 
cision) the numerical values are exactly the same. We 
suspect that ( Xbb/s) cfo = ( Xbb/s) pce but were un¬ 
able to show this analytically. It is also interesting to 
note that if we switch from the HRG model to the lattice 
data at around 160 MeV, then xbbT/s is approximately 
constant throughout the entire range 80 MeV <T< 220 
MeV, which is simple to deal with analytically. 

The BB susceptibility is also phenomenologically in¬ 
teresting, as it is important for net proton fluctuations 
which are used to search for the critical point. It is quite 
tempting to determine whether the inclusion of chemical 
freeze-out can improve the agreement of the HRG with 
the STAR measurements of these fluctuations [24, 25]. 
Simple considerations show that for any additional chem¬ 
ical potentials to have an effect on the skewness or kurto- 
sis, one must move beyond the Boltzmann approximation 
used in this paper. In our notation, the skewness mea¬ 
sured in experiment is 


Scr = T ®Xbb 
Xbb 


(43) 


Using (29), and noting that all baryons have Bf = 1, 


XbbT = ^ «i( T . Mi) ( 44 ) 

If = (45) 

' i 
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(a) (b) (c) 
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FIG. 3. (color online) Dimensionless susceptibilities in the case of full chemical equilibrium (FE), full freeze-out (FFO), and 
partial chemical equilibrium (PCE) assuming T c h = 160 MeV. The red points are lattice data adapted from [12]. 
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FIG. 4. (color online) Susceptibility to entropy density ratio in the case of full chemical equilibrium (FE), full freeze-out (FFO), 
and partial chemical equilibrium (PCE) assuming T c h = 160 MeV. The red points are lattice data adapted from [12]. 


Breaking the sum over i into baryons and anti-baryons 
and noting the only difference in the two is the sign of 
/is and Bi , we have 


9xbb 

dfi B 

XbbT 


— sinh(/r B /T) E ni(T,m) 

i, baryons 

2 cosh (^b/T) E n *( T ^) 

i, baryons 


Mb—0 


(46) 


0 


(47) 


And hence 


Sa = tanh(yUs/T) (48) 

which is independent of whether the system is in chemi¬ 
cal equilibrium or frozen out. Similarly, the kurtosis /ter 2 


can be shown to be equal to 1 whether or not one em¬ 
ploys chemical freeze-out. Thus, whether or not chemi¬ 
cal freeze-out of hadrons can improve the agreement of 
the HRG with STAR measurements requires more care¬ 
ful consideration which is beyond the scope of this paper; 
we defer it to future work. 


VI. CONCLUSION 


We have computed all six components of the suscepti¬ 
bility matrix for three conserved charges: electric charge, 























































































baryon number, and strangeness. We have done so using 
the hadron-resonance gas model employing two different 
implementations of chemical freeze-out. 

In almost all cases, the inclusion of chemical freeze-out 
enhances XafiT/s at low temperatures relative to the full 
equilibrium case. Thus we expect that the fluctuations 
at low temperatures will generally be larger than those 
found if one exclusively uses thermodynamic data from 
the lattice. This appears consistent with the findings 
of [25]. The exception is the diagonal electric charge 
susceptibility xqqT/s which is a decreasing function of 
T. Thus, the inclusion of chemical freeze-out will tend to 
reduce electric charge fluctuations relative to equilibrium. 

With regard to correlations and balance functions mea¬ 
sured in experiment [9, 10]: as compared to full equilib¬ 
rium, the inclusion of chemical freeze-out will tend to nar¬ 
row the pseudorapidity (Ar/) width of the balance func¬ 
tion (for a given diffusion constant). The reason is that 
enhanced fluctuations at late times (low temperatures) 
have less time to diffuse. The inclusion of resonance de¬ 
cays via PCE is unlikely to affect the balance functions 


except in the case of electric charge fluctuations. As seen 
in Fig. 4 (a), for xqqT/s , partial chemical equilibrium 
tends to significantly reduce the effect of chemical freeze- 
out, as the PCE curve lies closer to the full equilibrium 
curve than the one for full freeze-out. 

In summary, the inclusion of chemical freeze-out is an 
important difference between the thermodynamics of a 
physical heavy ion collision and lattice QCD calculations. 
This effect has implications for (among other things) par¬ 
ticle correlations and balance functions, and the search 
for the QCD critical point. Work to explore these impli¬ 
cations is already underway. 
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Appendix A: Particles, Resonances, and Decays 

We include all hadrons and resonances listed in the 
2014 PDG (with rating of *** or ****) with masses < 


2 GeV. The most massive particle we consider is the 
/2(2010), but we omit /o(600) (a meson) [26]. The 
masses and quantum numbers of all included particles 
match those given in the 2014 Review of Particle Physics 
provided by the PDG [22]. 

The measured decay modes and branching fractions of 
the higher mass resonances are very uncertain. Often 
decay channels are simply labeled as “seen” and/or the 
measured branching fractions do not sum to one. Hence, 
theoretical input is required. We rely almost exclusively 
on previous work by Josef Sollfrank and Pasi Huovinen 
who created a table of decays from the 2005 PDG data, 
supplemented by educated guesses based on the behav¬ 
ior of similar resonances and the fact that the branching 
ratios must sum to one [27, 28]. 

We have not attempted to update all branching frac¬ 
tions to incorporate more recent PDG data, as our results 
are most sensitive to lower mass hadrons/resonances, 
which were already measured precisely in 2005. Nev¬ 
ertheless, we have added three new resonances which 
were not present in the 2005 Review of Particle Physics: 
7T 2 (1880), A(1875) and JV(1900). The branching ratios 
for each of these resonances are very poorly known, hence 
(as in the original table) we rely on theoretical input and 
educated guesses to estimate decay rates of these parti¬ 
cles. These branching fractions do not significantly affect 
our results, but we include them here for completeness. 
Our assumptions are given in Table I, and we explain our 
rationale as follows. 

For the 7r 2 (1880), no decay modes are listed in the 
PDG. Our assumed branching fractions rely on the 
nearby resonance 7r 2 (1670) and the results of the model of 
[29]. For the A(1875), the PDG listings of the measured 
branching fractions do not sum to one. We attempt con¬ 
sistency with the measured values, forcing the sum to be 
one. For the _/V(1900), the measured branching fractions 
sum to less than 50%. We assume that the remainder is 
in the form of decays to Np , as is the case for the nearby 
resonance _/V(1720) (which has identical quantum num¬ 
bers). Finally, for the 04(2040), the PDG lists only a few 
decay modes as “seen”. Our input for this resonance is 
almost entirely theoretical, we start from the calculations 
of [30] and then include the KK and 'ifn modes to make 
the sum of the branching ratios one. 


tt 2 (1880) —► fm (33%) 

tt 2 (1880) —¥ / 2 (1270)tt (17%) 
7r 2 (1880) — ¥ pcu (17%) 

7r 2 (1880) —> KK * (11%) 

tt 2 (1880) — ¥ p( 1450)tt* (11%) 
tt 2 (1880) —¥ a(1320)?7 (11%) 
A(1875) —>~hhr (10 %) 

A(1875) —> Nlo (20 %) 

A(1875) —¥ Att (50 %) 

A(1875) —¥ Nttk (15 %) 
A(1875) —5- Np (5 %) 
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A(1900) - 

-> N TV 

(6%) 

A(1900) - 

-f pi1 

(12 %) 

IV(1900) - 

-»■ poj 

(12 %) 

A(1900) - 

AA 

(4%) 

A(1900) - 

->• E A 

(6%) 

A(1900) - 

->• Np 

(60 %) 

a 4 (2040) - 

—»• pu 

(37%) 

a 4 (2040) - 

—>■ p7T 

(20 %) 

a 4 (2040) - 

->• & 2 (1235)tt 

(22 %) 

o 4 (2040) - 

^/ 2 (1270)tt 

(11 %) 

o 4 (2040) - 

-> KK 

(5%) 

a 4 (2040) - 

—> rfn 

(5%) 


TABLE I: Assumed branching fractions for newly 
discovered resonances 


Appendix B: Sample Decay Coefficients for Partial 
Chemical Equilibrium 

As a convenience to the reader, Table II is a sample 
of the “effective number of stable particles” formed from 
each unstable particle (i.e. some of the dij coefficients). 
Note that this is only a partial list. 

While we have tabulated the dij coefficients for pho¬ 
tons produced in decays, we omit the photon when com¬ 
puting thermodynamic quantities like n, s and Xa/S as it 
is not expected that these photons will thermalize and 
they are not considered part of the HRG model. 
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7T + 

7T° 

7T 

A+ 

I<~ 

A 0 

A 0 

V 

UJ 

p 

p 

n 

n 

v' 

4> 

A 

A 

P + 

0 

1.000 

1.000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

P° 

0 

1.000 

0 

1.000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

P~ 

0 

0 

1.000 

1.000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

A'* (892)+ 

0 

0.667 

0.333 

0 

0.333 

0 

0.667 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

A* (892)” 

0 

0 

0.333 

0.667 

0 

0.333 

0 

0.667 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

A* (892)° 

0 

0 

0.333 

0.667 

0.667 

0 

0.333 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

A* (892)° 

0 

0.667 

0.333 

0 

0 

0.667 

0 

0.333 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

/o(980) 

0 

0.52 

0.52 

0.52 

0.11 

0.11 

0.11 

0.11 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a o (980)+ 

0 

0.844 

0 

0 

0.156 

0 

0 

0.156 

0.844 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a o (980)° 

0 

0 

0.844 

0 

0.078 

0.078 

0.078 

0.078 

0.844 

0 

0 

0 

0 

0 

0 

0 

0 

0 

cio (980) 

0 

0 

0 

0.844 

0 

0.156 

0.156 

0 

0.844 

0 

0 

0 

0 

0 

0 

0 

0 

0 

hi(1170) 

0 

1.000 

1.000 

1.000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

6i(1235) + 

0 

1.000 

0 

0 

0 

0 

0 

0 

0 

1.000 

0 

0 

0 

0 

0 

0 

0 

0 

&i(1235)° 

0 

0 

1.000 

0 

0 

0 

0 

0 

0 

1.000 

0 

0 

0 

0 

0 

0 

0 

0 

&i(1235)“ 

0 

0 

0 

1.000 

0 

0 

0 

0 

0 

1.000 

0 

0 

0 

0 

0 

0 

0 

0 

ai(1260) + 

0 

1.000.5 

1.000 

0.5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

ai(1260)° 

0 

1.000 

1.000 

1.000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

oi(1260)“ 

0 

0.5 

1.000 

1.5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

A++ 

0 

1.000 

0 

0 

0 

0 

0 

0 

0 

0 

1.000 

0 

0 

0 

0 

0 

0 

0 

A++ 

0 

0 

0 

1.000 

0 

0 

0 

0 

0 

0 

0 

1.000 

0 

0 

0 

0 

0 

0 

A+ 

0.006 

0.331 

0.663 

0 

0 

0 

0 

0 

0 

0 

0.669 

0 

0.331 

0 

0 

0 

0 

0 

A+ 

0.006 

0 

0.663 

0.331 

0 

0 

0 

0 

0 

0 

0 

0.669 

0 

0.331 

0 

0 

0 

0 

A 0 

0.006 

0 

0.663 

0.331 

0 

0 

0 

0 

0 

0 

0.331 

0 

0.669 

0 

0 

0 

0 

0 

A 0 

0.006 

0.331 

0.663 

0 

0 

0 

0 

0 

0 

0 

0 

0.331 

0 

0.669 

0 

0 

0 

0 

A“ 

0 

0 

0 

1.000 

0 

0 

0 

0 

0 

0 

0 

0 

1.000 

0 

0 

0 

0 

0 

A“ 

0 

1.000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1.000 

0 

0 

0 

0 


TABLE II: Incomplete list of the coefficients dij. The top row lists the 18 lightest stable particles. The first column lists the 
lightest 26 unstable particles. The table entry of the ith row, jth column is the coefficient dij. We omit columns for the heavier 
S, S, ft (their entries are all zero for the selection of unstable particles listed). 




















































